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Abstract 

(N 

U ' We consider the problem of partitioning the set of vertices of a given unit disk graph (UDG) 

D . into a minimum number of cliques. The problem is NP-hard and various constant factor ap- 

^1 ' proximations are known, with the current best ratio of 3. Our main result is a weakly robust 

polynomial time approximation scheme (PTAS) for UDGs expressed with edge-lengths, it either 
(i) computes a clique partition or (ii) gives a certificate that the graph is not a UDG; for the case 
(i) that it computes a clique partition, we show that it is guaranteed to be within (1 + e) ratio of 
\^ ' the optimum if the input is UDG; however if the input is not a UDG it either computes a clique 

partition as in case (i) with no guarantee on the quality of the clique partition or detects that it 
^ ' is not a UDG. Noting that recognition of UDG's is NP-hard even if we are given edge lengths, 

our PTAS is a weakly-robust algorithm. Our algorithm can be transformed into an O (^ujrfj 

time distributed PTAS. 

We consider a weighted version of the clique partition problem on vertex weighted UDGs 
i that generalizes the problem. We note some key distinctions with the unweighted version, where 

| ideas useful in obtaining a PTAS breakdown. Yet, surprisingly, it admits a (2+£)-approximation 

£SJ ■ algorithm for the weighted case where the graph is expressed, say, as an adjacency matrix. This 

s ! ' improves on the best known 8-approximation for the unweighted case for UDGs expressed in 

T^j- . standard form. 
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1 Introduction 



A standard network model for homogeneous networks is the unit disk graph (UDG). A graph 
G = (V, E) is a UDG if there is a mapping / : V ^ M 2 such that \\f(u) - f(v)\\ 2 < 1 ^ {u, v} G E; 
/(■ujj models the position of the node u while the unit disk centered at f(u) models the range of 
radio communication. Two nodes u and v are said to be able to directly communicate if they lie in 
the unit disks placed at each others' centers. There is a vast collection of literature on algorithmic 
problems studied on UDGs. See the survey [2]. 

Clustering of a set of points is an important subroutine in many algorithmic and practical 
applications and there are many kinds of clusterings depending upon the application. A typical 
objective in clustering is to minimize the number of "groups" such that each "group" (cluster) 
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1 /(.) is called a realization of G. Note that G may not come with a realization. 
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satisfies a set of criteria. Mutual proximity of points in a cluster is one such criterion, while points 
in a cluster forming a clique in the underlying network is an extreme form of mutual proximity. We 
study an optimization problem related to clustering, called the minimum clique partition problem 
on this UDGs. 

Minimum clique partition on unit disk graphs (MCP): Given a unit disk graph, G = 
(V,E), partition V into a smallest number of cliques. 

Despite being theoretically interesting, MCP has been useful for other problems. For example, 
[T7] shows how to use a small-sized clique partition of a UDG to construct a large collection of 
disjoint (almost) dominating sets. They [18] also show how to obtain a good quality realization 
of UDGs, and an important ingredient in their technique was to construct a small-sized clique 
partition of the graph. It is shown [12] how to use a small-sized clique partition to obtain sparse 
spanners with bounded dilation, which also permit guaranteed geographic routing on a related 
class of graphs. p3| employ MCP to obtain an 0(log* n) time distributed algorithm which is an 
0(logn)-approximation for the facility location problem on UDGs without geometry; they also 
give an 0(1) time distributed 0(l)-approximation to the facility location problem on UDGs with 
geometry also using MCP. Recently, |15j shows how to obtain a first O(l) approximation to the 
domatic partition problem on UDGs using MCP. 

On general graphs, the clique-partition problem is equivalent to the minimum graph coloring 
on the complement graph which is not approximable within n l ~ £ , for any e > 0, unless P=NP 
|22j . MCP has been studied for special graph classes. It is shown to be MaxSNP-hard for cubic 
graphs and NP-complete for planar cubic graphs [5]; they also give a 5/4-approximation algorithm 
for graphs with maximum degree at most 3. MCP is NP-hard for a subclass of UDGs, called unit 
coin graphs, where the interiors of the associated disks are pairwise disjoint [6]. Good approxi- 
mations, however, are possible on UDGs. The best known approximation is due to [6] who give 
a 3- approximation via a partitioning the vertices into co-comparability graphs, and solving the 
problem exactly on them. They give a 2-approximation algorithm for coin graphs. MCP has also 
been studied on UDGs expressed in standard form. For UDGs expressed in general form [18] give 
an 8-approximation algorithm. 

Our Results and Techniques: 

In this paper we present a weakly-robustd PTAS for MCP on a given UDG. For ease of exposition, 
first we prove this (in Section T2.ip when the UDG is given with a realization, /(.). The holy-grail is 
a PTAS when the UDG is expressed in standard form, say, as an adjacency matrix. However, falling 
short of proving this, we show (in Section [2]) how to get a PTAS when the input UDG is expressed 
in standard form along with associated edge-lengths corresponding to some (unknown) realization. 
The algorithm is weakly-robust in the sense that it either (i) computes a clique partition of the 
input graph or (ii) gives a certificate that the input graph is not a UDG. If the input is indeed a 
UDG then the algorithm returns a clique partition (case (i)) which is a (1 + ^-approximation (for a 
given e > 0). However, if the input is not a UDG, the algorithm either computes a clique partition 
but with no guarantee on the quality of the solution or returns that it is not a UDG. Therefore, this 

2 An algorithm is called robust if it either computes an answer or declares that the input is not from the restricted 
domain; if the algorithm computes an answer then it is correct [20]. We call our algorithm weakly-robust in that it 
always computes a clique partition or declares that the input is not from the restricted domain (i.e. not a UDG); if 
the input happens to be a UDG then the answer is a (1 + e)-approximate clique partition. Otherwise, it still returns 
a clique partition but there is no guarantee on the quality of the clique partition. 
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algorithm should be seen as a weakly-robust PTAS. The generation of a polynomial-sized certificate 
which proves why the input graph is not a UDG should be seen in the context of the negative result 
of [1] which says that even if edge lengths are given, UDG recognition is NP-hard. We show (in 
Section |4|) how this algorithm can be modified to run in O(^rr) distributed rounds. 

In Section [3] we explore a weighted version of MCP where we are given a vertex weighted UDG. 
In this formulation, the weight of a clique is the weight of a heaviest vertex in it, and the weight 
of a clique partition is the sum of the weights of the cliques in it. We note some key distinctions 
between the weighted and the unweighted versions of the problem and show that the ideas that 
help in obtaining a PTAS do not help in the weighted case. Yet, surprisingly, we show that the 
problem admits a (2 + ^-approximation algorithm for the weighted case using only adjacency. This 
result should be contrasted with the unweighted case where it is not clear as to how to remove the 
dependence on the use of edge-lengths, which was crucially exploited in deriving a PTAS. 

We use OPT to denote an optimum clique partition and opt to denote the size (or, in Section [3l 
weight) of an optimum clique partition. We also use n and m to denote the number of points (i.e. 
nodes of G = (V,E)) and the number of edges, respectively. 

2 A Weakly-Robust PTAS for UDG Expressed with Edge-lengths 

For simplicity, we first describe an algorithm when the input is given with a geometric realization. 
2.1 A PTAS for UDGs With a Geometric Realization 

We assume the input UDG is expressed with geometry of its points. Using a randomly shifted grid 
whose cell size is k x k (for k = k(e)) we partition the plane. Since the diameter of the convex hull 
of each clique is at most 1, for large values of A:, a fixed clique is cut by this grid (and therefore 
belongs to at most four cells) with probability at most |. Therefore, if we could efficiently compute 
an optimal clique partition in each k x k cell, then taking the union of these cliques yields a solution 
whose expected size at most (1 + e)opt. We can easily repeat this process O(logn) times to obtain 
a solution with size at most (1 + e)opt w.h.p. We call the algorithm MinCPl, formalized below. 

Theorem 1. Algorithm MinCPl (given below) returns, in poly-time, a clique partition of size at 
most (1 + e)opt w.h.p. 



Algorithm 1 MinCPl(G,e) 

1: Let k = ■ Place a grid whose squares have size k x k, on the plane. Call it Qofl- 
2: Pick (a, b) E [0, k) x [0, k) uniformly at random. 

3: Shift C?o,o by (<2) b) to get Q a ^ which is a grid shifted a units to the right and b units above. Q a ^ 

induces a random partition of V into points in k x k regions. 
4: for all k x k regions of Q a ^ do 

5: Obtain an optimal partition d for point-set P in the k x k square. 

6: Let C a: b = Ui=i Ci be the union of clique partitions obtained for the points in each k x k square. 
7: Repeat "Step 2-6" [logn] times and return the smallest C a ^ over the [logn] independent trials. 



We begin with a simple observation. 
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Observation 2. The diameter of the convex hull of every clique is at most 1. 



In the next subsection we argue how to perform "Step 5" of the algorithm MinCPl efficiently. 
Assuming this, we prove Theorem[TJ For a random shift Q a ^ and a clique C, we say that Q a ^ "cuts" 
C if some line of Q a t, crosses an edge of C. It is easy to see that: 



Pr 



C is cut 
by G a ,b 



< Pr 



a vertical or horizontal line of 5 a ,6 
crosses an edge of C 



2 

< - 

- k 



Thus, the expected number of cliques in an optimal partition that are "cut" by Q a \ 



2 -opt. So, by Markov's inequality, with probability at least 1/2 there are no more than 4 



is at most 

c UilCMl jr - Opt 

cliques cut by C/ a] b. Therefore, if we compute an optimal solution for each of the k x k grid cells and 
take the union of them, with probability at least 1/2 we get an excess of at most 4 x -r- opt cliques 
with respect to optimum since each clique that is "cut" by the grid can be counted up to four 
times. If we repeat this process for [logn] independent random trials, we get that with probability 
at least 1 — ^ the size of the solution we obtain is at most opt + opt < (1 + e)- opt. 



2.1.1 Optimal Clique Partition of a UDG in a k x k Square 

Unlike optimization problems such as maximum (weighted) independent set and minimum domi- 
nating set, where one can "guess" only a small-sized subset of points to obtain an optimal solution, 
the combinatorial complexity of any single clique in an optimal solution can be high. Therefore, it 
is unclear as to how to "guess" even few cliques, each of which may be large. A result of Capoyleas 
et al. [1] comes to our aid; a version of their result says that there exists an optimal clique partition 
where the convex hulls of the cliques are pair-wise non-overlapping. This phenomenon of separa- 
bility of an optimal partition, coupled with the fact that the size of an optimal partition in a small 
region is small, allows us to circumvent the above difficulty. The following simple lemma bounds 
the size of an optimal solution of an instance of bounded diameter. 

Lemma 3. Any set of points P in a k x k square has a clique partition of size 0(k 2 ). 

Proof. Place a grid whose cells have size 1/2 x 1/2. This grid induces a vertex partition where each 
block in the partition consists of the points that share a common grid cell (and therefore form a 
clique). □ 

We state a variant of a result by Capoyleas et al. [I] according to which there exists an optimal 
clique partition where the convex hulls of the cliques are non-overlapping, that is, for any pair of 
cliques in an optimal partition, there is a straight line which separates them| 

Theorem 4 Q4j). For a clique partition in which the convex hulls of the cliques are pairwise non- 
overlapping, there is a straight line lij that separates a pair of cliques Ci, Cj such that all vertices 
of Ci are on one side of lij, and all the vertices of Cj are on the other side of lij. (see Figure^ . 
Furthermore, this partition can be computed in poly-time. 

The general structure of the algorithm for computing optimal solution of a k x k cell is as 
follows. In order to reduce the search space for separator lines, one can find a characterization of 
the separator lines with some extra properties. Let Cj, Cj be a pair of cliques each having at least 

3 We gave a proof of this theorem [19] before it was brought to our attention that Capoyleas, Rote, and Woeginger 
[4] proved this much earlier in a different context. 
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(a) 



(b) 



Figure 1: (a) An optimal clique partition of UDG points in a bounded region; each light convex 
shape corresponds to a clique in the clique partition. The heavy line-segments represent segments 
of the corresponding separators, (b) A close-up view of Cj and Cj. A separator line, lij is shown 
which separates Cj and Cj, corresponding to the segment in (a). Note that l'^ is also a separator 
for Cj and Cj and 1'^ is passing through points x and y in Cj. 

two points. Let Lij be the (infinite) set of distinct separator lines. Since Ci and Cj are convex, 
there exists at least one line in that goes through two points of C, (or Cj) (see Figure Efb)). 
Therefore, given two cliques Ci and Cj in a clique partition (with pairwise non-overlapping parts) 
there is a separator line lij that goes through two vertices of one of them, say u,v G Ci such that 
all the vertices of Cj are on one side of this line and all the vertices of Ci are on the other side or on 
the line. Since there are 0{k 2 ) cliques in an optimal partition of k x k cell, there are 0(k 4 ) pairs of 
cliques in the partition and their convex hulls are pairwise non-overlapping. In fact, a more careful 
analysis shows that the dual graph of the regions is planar (see Figure Ufa)); thus there are 0(k 2 ) 
distinct straight lines, each of which separate a pair of cliques in our optimal solution. For every 
clique Ci, the separator lines (for all values of j) define a convex region that contains clique 
Ci. So once we guess this set of 0(k 2 ) lines, these convex regions define the cliques. We will try 
all possible (non-equivalent) sets of 0(k 2 ) separator lines and check if each of the convex regions 
indeed defines a clique and if we obtain a clique partition. This can be performed in 0(n k ) time 
(see [4] for more details). 



We weaken our assumption on having access to geometry; we assume only edge-lengths are known 
with respect to a feasible (unknown) realization of the UDG. We prove that, 

Theorem 5. Given a graph G with associated (rational) edge-lengths and e > 0, there is a poly- 
nomial time algorithm which either computes a clique partition of G or gives a certificate that G 
is not a UDG. If G is a UDG, the size of the clique partition computed is a (1 + e)- approximation 
of the optimum clique partition (but there is no guarantee on the size of the clique partition if the 
input graph is not UDG). 

The high level idea of the algorithm is as follows. As in the geometric case, we first decompose 
the graph into bounded diameter regions and show that if we can compute the optimum clique 
partition of each region then the union of these clique partitions is within (1 + e) fraction of the 
optimum. There are two main difficulties here for which we need new ideas. The first major 



2.2 A PTAS for UDGs With Edge-Lengths Only 
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difference is that, we cannot use the random shift argument as in the geometric case. To overcome 
this, we use a ball growing technique that yields bounded diameter regions. This is inspired by 
|13j who give local PTAS for weighted independent set, and minimum dominating set for UDGs 
without geometry. The second major difference is that, even if we have the set of points belonging 
to a bounded region (a ball) it is unclear as to how to use the separation theorem to obtain an 
optimal solution for this instance. Note that we are not guaranteed to have a UDG as input. We 
show that we can either compute a clique partition for each subgraph induced by a ball, or give a 
certificate that the subgraph is not UDG. If it is a UDG, then our clique partition is optimal but 
if it is not a UDG there is no guarantee on its size. 

Let B r (v) = {u : d(u, v) < r}, where by d(u, v) we mean the number of edges on a shortest path 
from u to v. So, B r (v) can be computed using a breadth-first search (BFS) tree rooted at v. We 
describe our decomposition algorithm which partitions the graph into bounded diameter subgraphs 
in Algorithm 2. We will describe a procedure, called OPT-CP which, given a graph induced by 
the vertices of B r (v) and a parameter I = poly(r), runs in time \B r (v)\°^ 2 ^ < and either 

produces a certificate that B r {v) is not a UDG or computes a clique partition of B r (v); this clique 
partition is optimum if B r (v) is a UDG. We only call this procedure for "small" values of r. 

Algorithm 2 MinCP2(G, e) 
1: C^0;/?^ [ C0 Il O gil;^ Cl /3 2 . 

{where cq is the constant in Lemma [9j and c\ is the constant in inequality ([I]).} 
2: while V / do 
3: Pick an arbitrary vertex v G V 
4: r <— 

{Let C r (v) denote a clique partition of B r (v) computed by calling OPT-CP} 
5: while \C r+2 {y)\ > (1 + £)■ \C r (v)\ do 
6: r <— r + 1 

7: if (r > p) or (OPT-CP {B r (v)) returns "not a UDG") then 

8: return "G is not a UDG" and produce B r (v) as the certificate 

9: e^C\JC r+2 (v) 

10: V<-V\B r+2 (v) 

11: return C as our clique partition 



Clearly, if the algorithm returns C on "Step 11", it is a clique partition. Let us assume that 
each ball B r (y) we consider induces a UDG and that the procedure OPT-CP returns an optimal 
clique partition C r (v) for ball B r (v). We show that in this case \C\ < (1 + e)opt. We also show that 
for any iteration of the outer "while-loop" , "Step 5" of MinCP2 is executed in time polynomial in 
n, by using edge-lengths instead of Euclidean coordinates. 

For an iteration i of the outer loop, let vi be the vertex chosen in "Step 3" and let r* be the 
value of r for which the "while-loop" on "Step 5" terminates, that is, \C r * + 2{vi)\ < ■ \C r * (vi)\. 

Let k be the maximum number of iterations of the outer loop. The following lemmas show that two 
distinct balls grown around vertices are far from each other, that the union of the optimal solutions 
to the balls form a lower-bound on the cost of the entire instance, and that the cost of C and opt 
is within a factor (1 + e) of opt. 

Lemma 6. For every i ^ j, every pair v E B r *(vi) and u £ B r *(vj) are non-adjacent. 
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Proof. Without loss of generality, let i < j. Therefore, every vertex in B r *(vj) is at a level larger 
than r* + 2 of the BFS tree rooted at vi, otherwise it would have been part of the ball B r * + 2[vi) 
thus removed from V . Note that in a BFS tree rooted at Vi, there cannot be an edge between a 
level r and r' with r 1 > r + 2. Thus there cannot be an edge between a node in v € B r *(vi), which 
has level at most r* and a node u 6 B r *(vj), which would been at a level at least r* + 3 in the BFS 
tree rooted at V{. □ 

Next, we derive a lower-bound on opt. 

k 

Lemma 7. opt > |C r *(t>j)| 
i=i 

Proof. Note that B r *(vi) is obtained by constructing a BFS tree rooted at vertex up to some 
depth r*. According to Lemma El there is no edge between any two nodes v £ B r *(vi) and 
u £ B r *(vj). So, no single clique in an optimum solution can contain vertices from distinct B r *(vi) 
and B r *{vj). Consider the subset of cliques in an optimal clique partition of G that intersect B r *(vi) 
and cafl this subset OPTj. The argument above shows that OPTj is disjoint from OPTj. Also, 
each OPTj contains all the vertices in B r *(vi). Since C r *(t>j) is an optimal clique partition for 
B r *(vi), |OPTj| > \C r *(vi)\. The lemma immediately follows by observing that OPTj and OPTj 
are disjoint. □ 

The next lemma relates the cost of our solution to opt. 

k 



Lemma 8. If \C r * + 2(vi)\ < (1 +e)- \C r *(vi)\, then 



\JC rt+2 (r 



i=l 



< (1 +£)• opt 



Proof. 



i=l 



^|C rr+2 (^)| < (l + eyJ2\C r *(vi)\ < (l+e)-opt □ 



i=l 



Finally, we show that the inner "while-loop" terminates in O(^), so r* £ O(j). Obviously, the 
"while-loop" on "Step 5" terminates eventually, so r* exists. By definition of r*, for all smaller 
values of r < r*: \C r (vi)\ > (1 +£)• \C r -2(vi)\. Since diameter of B r (vi) is O(r), if B r (v) is a UDG, 
there is a realization of it in which all the points fit into a r x r grid. Thus, |C r (uj)| G 0(r 2 ). So 
for some a S 0(1): 

a-r 2 > \C r ( Vi )\ > (l + e)-|C r _ 2 (t>i)| > ... > (1 + e)5- \C ( Vi )\ = 0({Vl + ~e) r ), 

when r is even (for odd values of r we obtain |C r (uj)| > (1 + e)~ ■ \C\{vi)\ > 0((y/l + e) r 1 ). 
Therefore we have: 

Lemma 9. There is a constant cq > such that for each i: r* < cq/e- log 1/e. 

In the next subsection, we show that the algorithm OPT-CP, given B r (v) and an upper bound 
£ on |C r (t>)|, either computes a clique partition or declares that the graph is not UDG; the size 
of the partition is optimal if B r (v) is a UDG. The algorithm runs in time n ^ 2 \ By the above 
arguments, if B r {v) is a UDG then, there is a constant c\ > such that: 

\C r (v)\ = 0(rf)< Cl - C f 2 W- e (1) 
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We can set i = |~ci^§ log 2 |] for any invocation of OPT-CP as an upper bound, where c\ is the 
constant in 0(r* 2 ). So, the running time of the algorithm is n olyl l ei \ 

2.3 An Optimal Clique Partition for B r {v) 

Here we present the algorithm OPT-CP that given B r {v) (henceforth referred to as G') and an 
upper bound I on the size of an optimal solution for G' , either computes a clique partition of it or 
detects that it is not a UDG; if G' is a UDG then the partition is optimal. The algorithm runs in 
time n°^\ Since, by Lemma [U i is a constant in each call to this algorithm, the running time of 
OPT-CP is polynomial in n. Our algorithm is based on the separation theorem [3]. Even though 
we do not have a realization of the nodes on the plane, assuming that G' is a UDG, we show how 
to apply the separation theorem [1] as in the geometric setting. We use node/point to refer to a 
vertex of G' and/or its corresponding point on the plane for some realization of G' . We will use 
the following technical lemma. 

Lemma 10. Suppose we have four mutually adjacent nodes p, a, b, r and their pairwise distances 
with respect to some realization on the Euclidean plane. Then there is a poly-time procedure that 
can decide if p and r are on the same side of the line that goes through a and b or are on different 
sides. 

Proof. First, we describe how to detect if the quadrilateral on these four points is convex or concave. 
If the quadrilateral is concave, then one of the points will be inside the triangle formed by the other 
three. There are three possible cases: r is inside, p is inside, or one of a or b is inside (see 
Figure [Hc)-(e)). There are four triangles each of which is over three of these four points. The 
quadrilateral is concave if the sum of the areas of three of these triangles is equal to the area of 
the fourth triangle. Equivalently, it is convex if sum of areas of two of the triangles is equal to the 
sum of areas of the other two. Given a triangle with edge lengths x, y, z, using Heron's formula, 
the area of the triangle is equal to y/2(x 2 y 2 + y 2 z 2 + z 2 x 2 ) — (x 4 + y 4 + z 4 ))/4. So the area of a 
triangle is of the form y/~A where A is a polynomial in terms of lengths of the edges of the triangle. 
Suppose that the areas of the four triangles over these four points are y/A~i, y/A?, y/A^, and \fA~i. 
We need to check if the sum of two is equal to the sum of the other two and we would like to 
do this without computing the square roots of numbers. For instance, suppose we want to verify 
y/Al+y/M = y/A^+x/A^. For this to hold, we must have A i + A 2 + 2y/A 1 A 2 = A 3 + A 4 + 2xJ A^A 4 . 
Verifying this is equivalent to verifying D + yj A\A 2 = y/ A3A4 where D = \(A\ + A 2 — A% — A4). 
Taking the square of both sides, we need to have D 2 + A\A 2 + 2DyJ A\A 2 = A3A4, which is the 
same as ^(^3^4 — D 2 — AiA 2 ) 2 = A\A 2 D 2 . Thus by comparing two polynomials of edge-lengths 
(and without computing square roots) we can check if the quadrilateral is convex or concave. 

Suppose the quadrilateral is convex. If r and p are on two opposite corners (see Figure [2ja)), 
then r and p are on different sides. In this case \rp\ + \ab\ > \ra\ + \bp\ and \rp\ + \ab\ > \rb\ + \ap\. 
If rp is one of its sides (see Figure E(b)), then \rp\ + \ab\ is not the largest of the above three pairs 
of sums. 

Now suppose that the quadrilateral is concave. The only case in which r and p are on two sides 
of line ab is when one of a or b is inside the triangle obtained by the other three (see Figure (2)^e) ) . 
In this case, the area of the largest triangle is the one that does not contain a or b. Thus, if we 
compute the square of the areas of the four triangle, we can detect this case too. □ 
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(a) (b) (c) (d) (e) 

Figure 2: The five non-isomorphic configurations needed to consider for a quadrilateral on four 
points in Lemma [TOl 

Assume that G is a UDG and has an optimum clique partition of size a < £. The cliques fall 
in two categories: small (having at most 2a — 2 points), and large (having at least 2a — 1 points). 
We focus only on finding the large cliques since it is easy to guess all the small cliques. Suppose 
for each pair Cj,C,- € OPT of large cliques, we guess their respective representatives, c, and Cj. 
Further, suppose that we also guess a separating line lij correctly which goes through points Uij 
and Vij. For a point p that is adjacent to c% or Cj we want to efficiently test if p is on the the same 
side of line ly as q (the positive side), or on Cj's side (the negative side), using only edge-lengths. 
Without loss of generality, let both and Vy belong to clique Cj. For every node p different from 
the representatives: 

• Suppose p is adjacent to all of Ci,Uij,Vij,Cj. Observe that we also have the edg and 
CiVij. Given the edge-lengths of all the six edges among the four vertices Ci,Uij,Vij,p using 
Lemma [TOl we can decide if in a realization of these four points, the line going through Uij } v%j 
separates the two points p and Cj or not. If p and Cj are on the same side, we say p is on the 
positive side of lij for Cj. Else, it is on the positive side of lij for Cj. 

• Suppose p is adjacent to Cj (and also to and Vij) but not to Cj. Given the edge-lengths of 
all the six edges among the four vertices c, , , v^j , p using Lemma [10] we can decide if in a 
realization of these four points, the line going through Uij,Vij separates the two points p and 
q or not. If p and Cj are on the same side, we say p is on the positive side of lij for Cj. Else, 
it is on the positive side of lij for Cj. 

For each Cj and all the lines lij, consider the set of nodes that are on the positive side of all 
these lines with respect to Cj; we place these nodes in Cj. After obtaining the large and the small 
cliques, we obtain sets C±, . . . , C a . At the end we check if each Cj forms a clique and if their union 
covers all the points. The number of guesses for representatives is and the number of guesses 

for the separator lines is n°^ a \ So there are a total of n°^ a ' configurations that we consider. 

Clearly, if G' is a UDG then some set of guesses is a correct one, allowing us to obtain an 
optimum clique partition. If G is not a UDG, we may still find a clique partition of G . However, 
if we fail to obtain a clique partition in our search then the subgraph is a certificate that G is not 
a UDG. 
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3 (2 + ^-Approximation for Weighted Clique Partition using Ad- 
jacency 



In this section we consider a generalization of the minimum clique partition on UDGs, which we call 
minimum weighted clique partition (MWCP). Given a node-weighted graph G(V,E) with vertex 
weight wt (v), the weight of a clique C is defined as the weight of the heaviest vertex in it. For a 
clique partition C = {C±, C2, ■ ■ ■ , C t }, the weight of C is defined as sum of the weights of the cliques 
in C, i.e. wt (C) = wt (U* =1 Ci) = Yl\=i w ^ iPi)- The problem is, given G in standard form, say, 
as an adjacency matrix, construct a clique partition C = {C±, C2, . . . , Ct} while minimizing wt (C). 
The weighted version of the problem as it is defined above has also been studied in different contexts. 
See [3 O [8] for study of weighted clique-partition on interval graphs and circular arc graphs. 

Observe that MWCP distinguishes itself from MCP in two important ways: (i) The separability 
property which was crucially used earlier to devise a PTAS does not hold in the weighted case, and 
(ii) the number of cliques in an optimal solution for a UDG in a region of bounded radius is not 
bounded by the diameter of the region anymore, i.e. it is easy to construct examples of weighted 
UDGs in a bounded region where an optimal weighted clique partition contains an unbounded 
(in terms of region diameter) number of cliques. In addition, examples where two cliques in an 
optimal solution are not separable, that is, their convex hulls overlap, is easy to construct. (See 
the examples given in Figure [3j) To the best of our knowledge, MWCP has not been investigated 
before on UDGs. We, however, note that a simple modification to the algorithm by [18] also yields 
a factor-8 approximation to the weighted case, a generalization which they do not consider. Here, 
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Figure 3: (a) Two overlapping weighted cliques, A = {a\, . . . , a^} and B = {b\, . . . , bk} are shown, 
ai,bi are independent for all i. The heavy polygon has vertices weighted k while the dashed ones 
are weighted 1. opt = k + 1 while any separable partition must pay a cost of at least 2k. (b) A 
UDG which is a matching between two cliques for which OPT contains t cliques. The weight is 
less than 2- a. 

we give an algorithm which runs in time 0(n poly ^ 1 / £ -') for a given e > and computes a (2 + e)- 
approximation to MWCP for UDGs expressed in standard form, for example, as an adjacency 
matrix. Our algorithm is weakly robust in that it either produces a clique partition or produces a 
polynomial-sized certificate proving that the input is not a UDG. When the input is a UDG, the 
algorithm returns a clique partition and it is guaranteed to be a (2 + e)-approximation; but if the 
input is not UDG there is no guarantee on the quality of the clique partition (if it computes one). 



Theorem 11. Given a graph G expressed in standard form, and e > 0, there is a polynomial 
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time algorithm which either computes a clique partition of G or gives a certificate that G is not a 
UDG. If G is a UDG, the weight of the clique partition computed is a (2 + e)- approximation of the 
minimum weighted clique partition (but there is no guarantee on the weight of the clique partition 
if the input graph is not UDG). 

Our algorithm will borrow some ideas developed in Section [2] and in [18]. The high level idea of 
the algorithm is as follows. Similar to the algorithm in Section [21 we first decompose the graph into 
bounded diameter regions and show that if we can compute a (2 + e)-approximate clique partition 
of each region then the union of these clique partitions is within (2 + e) fraction of opt. We will 
employ a similar ball growing technique (as in Section[2]) that will give us bounded diameter regions. 
We then show that we can either compute a clique partition or give a certificate that the subgraph 
is not a UDG. If the subgraph is a UDG, then our clique partition is within a factor (2 + e) of 
the optimal. For the case of bounded diameter region, although the optimum solution may have a 
large number of cliques, we can show that there is a clique partition with small number of cliques 
whose cost is within (1 + e)-factor of the optimum solution. First we describe the main algorithm. 
Then in Subsection 13.11 we show that for each subgraph B r (v) (of bounded diameter) there is a 
near optimal clique partition with 0(r 2 ) cliques. Then in Subsection 13 . 21 we show how to find such 
a near optimal clique partition. 

Let us denote the weight of the optimum clique partition of G by opt. As before, let B r {v) = 
{u : d(u,v) < r}, called the ball of (unweighted) distance r around v, be the set of vertices that 
are at most r hops from v in G. Our decomposition algorithm described below (see Algorithm 3) 
is similar to Algorithm 2 and partitions the graph into bounded diameter subgraphs below. The 
procedure CP, given a graph induced by the vertices of B r (v) and a parameter t = poly(r), runs in 
time n ^ 2 - 1 ) and either gives a certificate that B r (v) is not a UDG or computes a clique partition 
of B r (v); this clique partition is within a factor (2 + e) of the optimum if B r (v) is a UDG. We only 
call this procedure for constant values of r. In the following, let < 7 < V9+4E-3 ^ e a ra tional 
number. See Algorithm 3. 

Algorithm 3 MinCP(G,7) 

1: C^0;/3^ [ C0 Il O gIl;^ Cl /3 2 . 

{where cq is the constant in Lemma [Tol and c\ is the constant in inequality ©.} 
2: while V / do 
3: v <— arg max. u {wt (u)} 
4: r <— 

{Let C r (v) denote a factor- (2 + 7) partition of B r (v) computed by calling CP} 
5: while wt (C r+2 (v)) > (1 + 7)- wt (C r (v)) do 
6: r <— r + 1 

7: if (r > p) or (CP{B r (v),£) returns "not a UDG") then 

8: return U G is not a UDG" and produce B r (v) as the certificate 

9: C^C\JC r+2 (v) 

10: V^V\B r+2 (v) 

11: return C as our clique partition 



Let k is the maximum number of iterations of the outer "while-loop" . The proof of the following 
Lemma is identical to the proof of Lemma [6j 

Lemma 12. Every two vertices v £ B r *(vi) and u G B r *(vj) are non-adjacent. 
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The following lemma shows a lower-bound for opt. 
Lemma 13. (2 + 7)- opt > wt Uj C r »(«i)^ 

Proof. Note that B r *(vi) is obtained by constructing a BFS tree rooted at vertex Vi up to some 
depth r*. Since the algorithm removes a super-set, B r * + 2(vi), which has two more levels of the BFS 
tree, using the previous lemma there is no edge between any two nodes v S B r *(vi) and u S B r *(vj) 
for any pair i 7^ j. So, no single clique in an optimum solution can contain vertices from distinct 
B r * (yi) and B r * (vj). Consider the subset of cliques in an optimal clique partition of G that intersect 
B r *(vi) and call this subset OPTj. The argument above shows that OPTj is disjoint from OPT.,. 
Also, each OPTj contains all the vertices in B r *(vi). Since C r *(vi) is a factor- (2 + 7) approximation 
for B r * (vi), (2 + 7)- wt (OPTj) > wt (C r * • The lemma immediately follows by observing that 
OPTj' and OPTj are disjoint. □ 

We can relate the cost of our clique partition to opt as follows. 
Lemma 14. If wt (C r * +2(^1)) < (1 + r y)wt(C r * (i>i)), then wt {\^ CV*+2(^i)^ < (2 + e)opt. 




= ]Twt (CV* +2 «)) < (l+ 7 )-^wt(C r .(0) < (2 + 7)(l + 7 )-opt, 

i=l i=l 

where the last inequality uses Lemma [131 D 

Next, we show that the inner "while-loop" terminates in O(-), that is each r* is bounded 

by O(^). This is similar to the proof of Lemma [9l Since the while loop terminates, r* exists 
and by definition of r* , it must be the case that for all smaller values of r < r*, wt (C r (t>j)) > 
(1 + 7)- wt (C r _2(i>i)). Because the diameter of B r (vi) is 0(r), if B r (v) is a UDG, there is a 
realization of it in which all the points fit into a r x r grid. Also, since t> j is a heaviest vertex in 
the (residual) graph, there is a clique partition whose weight is at most a- wt (uj) -r 2 . Therefore, 
wt (C r {vi)) < a- wt (vi) ■ r 2 , for some constant a. So: 

a-wt (uj)-r 2 > wt(CVK)) > (1+7)- wt (C r -2K)) > ••• > (1+7)^ wt (Cofa)) = wt • ( V 1 + tV , 



which implies a-r 2 > (^1 + jY , for the case that r is even. If r is odd we obtain a-r 2 > 
(■\/l + 7)^ 1 . Thus, the following lemma easily follows: 

Lemma 15. There is a constant Co > such that for each i: r* < 00/7- log I/7. 

In Subsection [321 we show the algorithm CP that given B r (v) and an upper bound £ on |C r (v)|, 
either computes a clique partition (which is within a factor 2 + 7 of opt if B r (v) is a UDG) or detects 
that the graph is not UDG; the algorithm runs in time >. By the above arguments, if B r (v) is 
a unit disk graph then there is a constant c\ > such that: 

|a( V )| = 0(r* 2 )< Cl -4log 2 - (2) 

T 7 
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We can set I = \c\-% log 2 ~] for any invocation of OPT-CP as an upper bound, where c\ is the 
constant in 0(r* 2 ). So, the running time of the algorithm is n ^ 1 ^ 4 ^. 

3.1 Existence of a Small Clique Partition of B r (v) having Near-optimal Weight 

Unlike the unweighted case, an optimal weighted clique partition in a small region may contain a 
large number of cliques. Yet, there exists a partition whose weight is within a factor (1 + 5r) of the 
minimum weight which contains few cliques (where by "few" we mean i as in Algorithm 3). The 
existence of a light and small partition allows us to enumerate them in the same manner in the 
algorithm of subsection 12.31 yielding a (2 + 7)-approximation for the problem instance in a ball of 
small radius. In the following, let r £ O(^); we focus on the subproblem that lies in some B r (v). 
Recall that any ball of radius r can be partitioned into 0(r 2 ) cliques (Lemma [3j). We begin with a 
simple lemma which states that for any clique partition C, if the set of vertices can be be covered 
by another clique partition C containing x cliques then the sum of the weights of the x cliques in 
C' is not significantly more than the weight of the heaviest clique in C. 

Lemma 16. For any collection of disjoint cliques C = {Ci,C2, • • • ,Ct} having weights such that 
wt(C\) > wt(C2) > ■ • • > wt{Ct) suppose the vertices of C can be partitioned into x cliques 
C' = {C[,C' 2 ,...,C X }. Then wt(C) = wt(\J l=l C[) = Ei=i wt ( c i) < x-wt(d) 

Proof. Without loss of generality, let wt(CJ) > wt(C' 2 ) > ... > wt(C' x ). Since C partitions 
vertices in C, wt (C[) = wt (Ci). Since \C'\ = x, wt (C) = £f =1 wt (C[) < x- wt (C[) = x- wt (Ci). 

□ 

In an optimal partition of a ball of radius r, the sum of the weights of the lighter cliques is not 
significantly more than its weight. 

Lemma 17. Let C = {C±, C2, . . . , C{\ be an optimal clique partition and let wt(Ci) > wt (C2) > 
... > wt(Ct). Suppose there is another clique partition C = {C[, . . . ,C' X } of the vertices of C. 
Then, for every 1 < i < t: (x - 1)- wt(d) > E*=i+i wt ( c i)- 

Proof. By way of contradiction, suppose there exists an index 1 < j < t such that (x — 1)- wt (Cj) < 
Yli=j+i w * (Ci)- Because |Jz=i C can be covered by C, so can U/=j C Let the 2 < x' < x be 
the smallest index such that Cj = {C[ , C 2 , ■ ■ ■ , C' x ,} covers |j| =J Ci. On the other hand, (x' — 
1)- wt (Cj) <(x- 1)- wt (Cj) < E*=j+i wt (C), which implies 

< j j'-l 

opt = wt (Ci) > wt + ( x ' ~ X )' wt ( C J') = wt (C) + x'wt (Cj) . (3) 
1=1 1=1 1=1 

By LemmafTBl wt [C'A < x' • wt (Cj). This, combined with inequality ([3]) implies opt > Yli=i w ^ (C) + 

Ef=i w * Therefore the cliques in C" = {C±, C2, ■ ■ ■ , Cj-\, C[, C 2 , • • • , C' x ,} cover all the nodes 
of cliques in C and has cost smaller than opt. If a vertex belongs to two or more cliques in C" we 
remove it from all but one of them to obtain a clique partition with cost no more than cost of C" 
which is smaller than opt. This completes the proof. □ 
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We now are ready to prove the main result of this section which states that for any optimal 
weighted clique partition of a ball of radius r, there exists another clique partition whose weight 
is arbitrarily close to the weight of the optimal partition, but has 0(r 2 ) cliques in it. Since the 
radius of the ball within which the subproblem lies is small, r G O(-), this means that if we were 
to enumerate all the clique partitions of the subproblem up to 0(r 2 ), we will see one whose weight 
is arbitrarily close to the weight of an optimal clique. Choosing a lightest one from amongst all 
such cliques guarantees that we will choose a one whose weight is arbitrarily close to the optimal 
weight. 

Lemma 18. Let 7 > and r G 0(1/7) be two constants. Let C = {C\, O2, . . . , Ct} be an optimal 
weighted clique partition of B r (v) and let C = {C[, . . . , C' x } be another clique partition of vertices 
of C with x G 0(r 2 ). Let wt(C\) > wt(C2) > ■ ■ ■ > wt(Ct). Then, there is a partition of vertices 
of C into at most j + x cliques for some constant j = with cost at most (1 + ^)opt. 

Proof. Without loss of generality, we assume that both x and t are at least two (as if B r (v) is a 
clique we are done). Consider an arbitrary value of j < t. Since Uz=i Ci can be covered by x cliques 
in C 1 , there is an index x' (2 < x' < x) such that U*=j &i can be covered by Cj = {C[, C 2 , . . . , C' x ,}. 
By applying Lemma [T71 repeatedly: 



opt > wt (Q) > wt (Q) ) +X^wt(C/) > ... > ( -^—) ' 1 • wt {C 3 



1=1 x 



( x ' _ 1 y'-i 

opt 1 > >a/-wt(Cj) (4) 



Using inequality (j3J): 

opt + ° Pt " {X ',~ 2 1)J l >E wt (^) + wt (CT,-) > J2 wt w + E wt ( c 



x'3- 

1=1 1=1 1=1 



(5) 



where the second inequality follows by applying Lemma [TBI Let C" = {Ci, . . . , Cj-i, C[, . . . , C' x ,}. 
Thus, the cliques in C" cover all the vertices of C and has total cost at most ( 1 + — ~,}-2 — ) °pt by 



inequality ([5|) . If a vertex belongs to two or more cliques in C" we remove it from all but one of them 
arbitrarily to obtain a clique partition of size j — 1 + x' and whose total cost is upper bounded by 

(l + ^l" 1 ) opt. Note that, ^''PS' = W ~ l ) (^ i )' ? ~ 2 and < ^ < 1 (because x' > 2). 
Since r G O(^) and x' < x G 0(r 2 ), for an appropriate choice of j = .7(7), (x' — 1) ( ^-7^- ) < 7/2. 



Thus we obtain a clique partition with j + x — 1 cliques and cost at most (1 + 7/2)- opt. This proves 
the lemma. □ 

3.2 (2 + 7)-Approximation for MWCP in B r (v) 

Finally, we show how to compute a (2 + 7)-approximate MWCP of the graph B r (v) for any given 
7. For an edge ordering L = (e±, e2, • • • , e m ) of a graph G with m edges, let Gl[i\ denote the 
edge induced subgraph with edge-set {ei,ej+i, . . . ,e m }. For each e,, let Nl[i] denote the common 
neighborhood of the end-points of ej in An edge ordering L = (ei, e2, . . . , e m ) is a CNEEO if 
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for every ej in L, Ni{i) induces a co-bipartite graph in G. It is known [20] that every UDG graph 
admits a co-bipartite edge elimination ordering (CNEEO). In the following, let G v denote B r {v). 
We state a lemma of [18]. 

Lemma 19. 118)1 Let C be a clique in G v , and let L be a CNEEO of G v . Then, there is an i, 
1 < i < m, such that Ni[i] contains C. 

Assume that G v can be partitioned into a < I = 0(l/7 2 ) cliques, O = {0±, O2, • • • , O a }, such 
thatwt(O) < (l + 2)-wt(OPT„), where OPT v is an optimal weighted clique partition of G v . Note 
that by Lemma [TBI this is true for subgraph B r {y). Suppose that we are given the upper bound I] we 
will try all possible values of a. Without loss of generality, let wt (0±) > wt (O2) > . . . > wt (O a ). 
Observe that, without loss of generality, we can assume Oi is a maximal clique in [JJ =i Oj. The 
implication of the above lemma is that even though we do not know Oi, hence we do not know 
0, we do know that for every CNEEO L of G v , there is an ej such that Nl[i\ can be partitioned 
into at most two cliques that fully cover 0\. Since 0\ is a heaviest clique, the two cliques that 
cover the subgraph Nl[i\ pay a cost of at most 2- wt (Oi). This suggests an algorithm that guesses 
an edge sequence /2, . . . , f a ) of G v . Then, the algorithm computes L, a CNEEO of G v . The 
algorithm's first guess is "good" if f\ is an edge in 0\ that occurs first in L. Suppose that this is 
the case and suppose that f\ has rank i in L. Then, 0\ is contained in iVx,[i], and we cover Nl[i\ 
with at most two cliques. Call these C[ and C" and wt (C[) + wt (G") < 2- wt (Oi). So, when we 
remove Nl[i\ from G v , we get a UDG which can be partitioned into at most a — 1 cliques, namely, 
O' = {0 2 , ■ ■ ■ , O a }. We then again construct a CNEEO, L' , of G' v = G v \ N L [i]. Just like before, 
our guess f% is "good" if fa is an edge in O2 and occurs first in L' . Let i 1 be the rank of /2 in L' , we 
see that Ni,i[i'] fully contains O2, and we again cover it with at most 2 cliques. Next, delete Nii\i'\ 
from G' v to get a graph which can be partitioned into a — 2 cliques, and so on. See Algorithm 4 
for details. 

Algorithm 4 CP(G V ,£) 
1: C <— V; min <— wt (C); 
2: for all a < £ do 

3: for all a-edge sequence (/1, /2, ■ • ■ , /«) of G v do 

4: Go *- G v 

5: for j = 1 to a do 

6: Compute a CNEEO L of Gj_i 

7: % <— rank of fj in L 

8: Partition Nl[i] into two cliques C'j and Cj' 

9: Gj <— Gj_i \ A r L[i] 

10: if G Q = and wt (\J" =1 {C' jf C'!}\ < min then 

11: C <- UJ=i{Cj, G;|; min <- wt (C); 

12: return C 



Note that while Lemma [19] allows us to cover any clique with at most 2 cliques, it does not 
find the clique. In the algorithm, note that if at any point, the algorithm is unable to construct 
a CNEEO, we can declare that the graph G v is not a UDG. Also, if for all invocations of the 
algorithm by an external algorithm that guesses the value of opt we are unable to find a clique 
partition, then again we can declare that G v is not a UDG. 
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4 0(log* n)-round Distributed PTAS for UDGs with Edge-Lengths 



In this section, we give details of a distributed PTAS for MCP which runs in 

0(%rr) rounds of 

distributed computation under the LOCAL model of computation [IB]. The model of computation 
that we employ assumes a synchronous system where communication between neighboring nodes 
takes place in synchronous rounds using messages of unbounded size [16]. So, in a single round 
of communication, any node acquires the subgraph (information pertaining to the set of nodes, 
edges, the states of local variables, etc.) within its immediate neighborhood. So, after k rounds of 
communication, any node acquires complete knowledge about its /c-neighborhood. 

Observe that in Algorithm MinCP2, the radius r of any ball B v (r) is bounded above by 0(l/e), 
while the center, v, is an arbitrary vertex. Since the radius of any ball is "small", the maximum 
number of rounds of distributed computation that the sequential algorithm needs before terminating 
the "while- loop" is also "small". Therefore, for any pair of balls B u (ri) and B v (rj), such that 
d(u,v) £ uj(l/e), one should be able to run part of the sequential algorithm in parallel, as they 
surely are independent of each other. We borrow some ideas from [TO] and find regions that are far 
apart such that we can run the sequential algorithm in those regions in parallel. See Algorithm 5 
for details. 

Algorithm 5 Distr-MCP-UDG(G, e) 

1: (3 <— [co^ log i] ; i <— ci/3 2 ; all vertices are unmarked. 

{co is the constant in Lemma [9] and c\ is the constant inequality ([1]).} 
2: Construct a maximal subset, V c C V, such that for any pair u, v G V c , d(u,v) > (3. Construct 

a graph G c = (V C ,E C ), where E c = {{u, v} : u,v £ V c , dc(u,v) < 4/?}. We call V c , the set of 

leaders. 

3: Proper color G c using A(G C ) + 1 colors, where A(G C ) is the maximum degree of G c . 

4: Every v 6 V \V C , "assigns" itself to a nearest leader u £ V c , with ties broken arbitrarily, and 

colors itself the same color as the leader. 
5: for i = 1 to A(G C ) + 1 do 

6: For each leader j with color i let G\ be the subgraph induced by the vertices assigned to 

leader j. 
7: for all G\ in parallel do 

8: Consider a fixed ordering on the unmarked vertices of G\ ; 

9: Run the sequential ball growing algorithm on the next (in this ordering) unmarked vertex 
f £ G], we compute B r (v); Note that B r {v) might contain vertices of different colors (from 
outside G\). 

10: Compute (using the sequential algorithm) the optimal clique-partition of B r (v) and "mark" 
all those vertices 



It should be pointed out that adapting the algorithm of [TO] for maximum independent set and 
minimum dominating set to our setting is not trivial. The reason is that MCP is a partition of the 
entire vertex set and partitioning just a subset well enough will not do. Specifically [TO] chooses a 
subset of vertices upon which their ball-growing algorithm is run; it suffices for their purposes to 
dispense with the remaining subset of vertices that were not picked by their ball-growing algorithm. 
If we had followed a similar scheme then we would surely get a good clique partition on a subset 
of vertices; however, it is unclear as to how to obtain a good partition of the remaining subset in 
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terms of the optimal size for the original problem instance over the entire vertex set. As a means 
to circumvent this issue, we first construct a "crude" partition of the vertex set, instead of just a 
subset of vertices as done in |10| . 

4.1 Analysis 

We now show that the algorithm constructs a (1 + ^-approximation to MCP on UDGs given only 
rational edge-lengths in 0(log* n) rounds of distributed computation under the LOCAL model; we 
first show correctness of the algorithm, followed by bounding the number of communication rounds. 

Correctness: We prove that our algorithm is correct by showing that any execution of Distr- 
MCP-UDG can be turned into a sequential execution of MinCP2. As stated, every vertex has the 
same color as its leader; let a leader vertex be its own leader. We first show that the distance of 
every vertex to its leader is small. 

Lemma 20. For any vertex v ^ V c , there is a vertex u G V c such that dc(u,v) < (3. 

Proof. Suppose not. So there is a v whose distance to every u £ V c is more than (3. But then 
V' c = V c U {v} has the property that for all x,y £ V' c , dc(x,y) > (3, contradicting the maximality 
of V c . □ 

Next we show that for any pair of vertices u, v of the same color but with different leaders, the 
minimum distance between them is large enough so that a pair of balls of radius at most (3 over 
them will be disjoint, where [3 is defined in MinCP2 and Distr-MCP-UDG. 

Lemma 21. Consider two leaders x,y of the same color, say i, and any two vertices u 6 G'f and 
v £ G\ (note that we might have u = x or v = y). Then for all values of r considered in the ball 
growing algorithm, B r {u) and B r (v) are disjoint. 

Proof. Since x and y have the same color d(x, y) > 4/3. By Lemma [20l any vertex in either of Gf or 
G\ is at a distance of at most (3 from the respective leader; so, dc(u,v) > 2(3. The lemma follows 
easily by noting the fact that r < (3 in the ball growing algorithm. □ 

We are now ready to prove the correctness of Distr-MCP-UDG by showing an equivalence 
between any execution of it to some execution of MinCP2. 

Lemma 22. Any execution of Distr-MCP-UDG from "Step 5" to "Step 10" can be converted to a 
valid execution of MinCP2. 

Proof. Consider an arbitrary execution of Distr-MCP-UDG. Suppose that Vi,V2,Vs, . . . is a se- 
quence of disjoint sets of the vertices of V such that we run the ball growing algorithm in parallel 
(during Distr-MCP-UDG) on vertices of V\ (and thus we compute an optimal clique partition on 
each vertex of V\ in parallel) then we do this for vertices in V2, and so on. Note that the vertices 
in Vi all have the same color and each has a different leader. Consider an arbitrary ordering 7Tj 
of the vertices in each Vi and suppose that we run MinCP2 algorithm on vertices of V\ based on 
ordering tti, then on vertices of Vi based on ordering TT2, and so on. Since the vertices in each Vi 
have distinct leaders, by Lemma [2T1 the balls grown around them are disjoint. It should be easy 
to see that the balls grown by algorithm MinCP2 is exactly the same as the ones computed by 
Distr-MCP-UDG. □ 
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The following result follows immediately as a corollary to Lemma [22l 

Corollary 23. Given an e > 0, Distr-MCP-UDG constructs a clique partition of the input graph 
G with associated edge-lengths, or produces a certificate that G is not a UDG. If G is a UDG then 
the size of the partition is within (1 + e) of the optimum clique partition. 

Running Time: We now show that the algorithm runs in O( 1 °o( 1 " ) distributed rounds under the 
LOCAL model of computation. 

Lemma 24. "Step 2" requires 0(/?-log*n) rounds of communication. 

Proof. Observe that the result of "Step 2" is identical to constructing a maximal independent set 
(MIS) in Note that G 13 is also a UDG where the new unit is (3. As a result, G 13 is a subclass 
of growth-bounded graphs [5] where all the distances are scaled by (3; computation of MIS on G 3 
takes 0((3- log* n) rounds [21] while the construction of G@ takes (3 rounds. Hence, the number of 
rounds needed by "Step 2" can be bounded by 0{(3- log* n). □ 

It is easy to see that constructing G c requires at most 4/3 communication rounds. Next, we 
show that the maximum degree of G c , A(G C ) is bounded by a constant. 

Lemma 25. A(G C ) G O(l) 

Proof. Let v be a vertex of G c having maximum degree. Note that all its neighboring vertices in G c 
lie in a disk of radius at most 4/3. Also note that due to "Step 2" the minimum distance between 
any pair of vertices in G c is more than (3. As a result, any disk of diameter (3 contains at most 1 
vertex of G c . Using standard packing arguments of the underlying space, a crude upper bound on 
the number of vertices of G c in a disk of radius at most 4/3 is 256 vertices; this also upper bounds 
the degree of v. □ 

Next, we bound the number of rounds needed for "Step 3" 

Lemma 26. "Step 3" requires 0(/3-log*n) rounds of communication. 

Proof. For graphs whose maximum degree is A, a A + 1 proper coloring requires 0(A + log* n) 
rounds [11]. Since A(G C ) G 0(1) (Lemma [25]) . and the fact that distances in G c are scaled by a 
factor of 4/3 as compared to the distances in G, a A(G C ) + 1 proper coloring of G c can be obtained 
in 0(/3-log*n) rounds. □ 

"Step 4" requires at most (3 rounds of communication; according to Lemma 1201 for every 
v ^ V c , there is some u G V c that is at a distance at most (3 from it. The identity and color 
of such a vertex can be obtained in (3 rounds. We can now bound the number of rounds that 
Distr-MCP-UDG requires. First, note that for any iteration, i, of "Step 7", only knowledge of a 
subgraph up to radius (3 is required, and any node can obtain knowledge of the subgraph up to 
radius (3 from it in (3 rounds of communication. So, for any vertex in G J - obtains knowledge about 
the "marked/unmarked" status of all the vertices in G\ in (3 rounds of communication. Since the 
diameter of each G\ is at most 2/3, the number of balls to grow in "Step 9." is at most 0((3 2 ). 
Therefore: 

Theorem 27. Distr-MCP-UDG requires 0(/3-log* n) rounds of communication under the LOCAL 
model of computation. 
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5 Concluding Remarks 



Recall that the weakest assumption that we needed to obtain a PTAS for unweighted clique partition 
problem was that all the edge lengths are given. This information was crucially used in obtaining 
a robust PTAS. In the case of weighted clique partition, we gave a (2 + ^-approximation algorithm 
without the use of edge-lengths (using only the adjacency information). It will be interesting to 
see if a PTAS exists for the unweighted case but with reliance only on adjacency. 

It is also unclear if a PTAS is possible even with the use of geometry in the weighted case. 
Recall that the PTAS given in Sections [2] crucially uses the idea of separability of an optimal clique 
partition. However, in the weighted case, even though a near optimal clique partition in a small 
region has few cliques, there are examples where any separable partition pays a cost at least factor-2 
to that of a near optimal partition. We give an example in Figure El^a). In the example shown 
in Figure [3ja) two cliques of optimal weight are shown: one of them, A, whose vertices are the 
vertices of the A;-gon shown in dashed-heavy lines, and the other, B, whose vertices are the vertices 
of the fc-gon shown in solid-heavy lines. The example is that for k = 7. The vertices of A are 
labeled 01,02, . . . in a counter-clockwise fashion. The vertices of B are labeled such that 6j is 
diametrically opposite to Oj. The distance between a% and hi is more than 1 while the distance 
between at and bj, i ^ j is at most 1. So, there is an edge between m to every ai and to every 
7^ i- This is also the case for b{. In the figure, the edges incident to a\ are shown by solid-light 
lines. Also, the dashed arc shows part of the unit disk boundary that is centered at a\ - note that 
it does not include b±. Let the weights of vertices in A be k and the weights of vertices in B be 1. 
Clearly, opt < k + 1. However, any separable clique partition pays a cost of at least 2k: if vertices 
in A must all belong to a common clique, then every vertex in B must belong to a distinct clique 
in a separable clique partition. Also, note that as-per separability, a line going through {pi,P2} 
separates two cliques having weight 2k also. 

Note that our results only apply in the Euclidean plane; they do not generalize. In particular, 
Capoyleas et al. [I] give an "unseparable" instance in M 3 . Our result in the weighted case also is 
restricted to the plane; the concept of co-bipartite neighborhood edge elimination ordering (CNEEO) 
does not generalize to M 3 . 
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